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Abstract. 


The  paper  analyzes  in  detail  the  problem  of  various  plate  models, 
properties  of  their  solution  and  the  question  how  well  these  models 
approximate  the  3-dimensional  formulation.  The  boundary  layer  and  corner 
singularities  of  the  solution  are  analyzed.  A  hierarchy  of  the  models 
converging  to  the  3-dimensional  solution  is  constructed.  The  detailed 
numerical  computations  show  various  basic  aspects  of  plate  modeling  in  a 
concrete  setting. 

Key  Words:  Plates,  Kirchhoff  model,  Reissner-Mindlin  model,  plate 

paradox,  corner  singularities 

1.  Introduction. 

The  plate  and  shells  are  basic  elements  in  structural  mechanics. 

Historically  much  attention  has  been  given  to  the  derivation  of  plate  and 
shell  models  which  lead  to  the  different  solution  depending  on  the  model 
used;  see  e.g.  [1]  for  a  survey.  The  derivation  of  a  model  is  typically 
based  on  various  mechanical  considerations  and  principles.  Only  in  1959  in 
[2]  was  the  first  rigorous  proof  of  the  relation  between  the 
three-dimensional  solution  and  the  plate  model  given. 

Many  papers  then  had  analyzed  the  plate  theories  via  asymptotic 
arrivals.  See  e.g.,  [3]  and  [4]  and  references  there.  Current  development 
of  numerical  methods  leads  to  the  understanding  of  the  plate  and  shell 
modeling  as  the  application  of  the  dimensional  reduction  principles  to  the 
three-dimensional  problem.  This  then  loads  to  a  hierarchical  modeling  and  to  “ 

a  convergent  sequence  of  models.  Adaptive  approaches  together  with  the  — - - 

a-posteriori  analysis  of  the  error  of  the  model  is  a  realistic  goal  to  be  '  Codes 


achieved. 


The  recent  development  of  the  h-p  version  of  the  finite  element  method 
allows  to  create  a  natural  hierarchy  of  models  based  on  the  polynomial 
approximation  through  the  thickness  and  the  use  of  the  variational 
approaches.  Such  approximation  is  in  a  certain  sense  optimal.  See  e.g.  [5]. 
The  h-p  version  also  allows  for  assessing  the  accuracy  of  the  finite  element 
solution  when  compared  with  the  exact  solution  of  the  plate  model  under 
consideration,  as  well  as  the  error  when  compared  with  the  three-dimensional 
solution. 

In  this  paper  we  discuss  the  accuracy  of  various  plate  models  when 
compared  with  the  solution  of  the  three-dimensional  problem.  Although  the 
3-dimensional  solution  is  not  exactly  known,  we  can  find  it  by  numerical 
means.  The  error  lies  in  a  tolerance  bracket  which  is  adjusted  to  the  aims 
of  our  analysis.  Hence  we  can  assume  that  the  3-dimensional  solution  is 
available.  The  accuracy  control  is  based  on  the  comparison  of  the  solution 
for  various  meshes  and  degrees  of  elements.  The  computation  was  made  by  the 
program  MSC/PROBE  with  h-p  capabilities. 

Typically  we  deal  with  a  uniformly  loaded  square  plate,  |x^|  <0.5, 
i  =  1,2,  with  the  plate  thickness  d  *  .  Because  of  symmetries,  the 

problem  is  solved  for  the  quarter  only.  The  basic  mesh  in  the  quarter  of  the 
plate  is  defined  by  the  lines  =  x2  ~  ^i  various  meshes 

1)  25  elements: 

^  =  0.  0.36,  0.492.  0.4994,  0.49988,  0.5 

2)  36  elements: 

^  =  0,  0.3,  0.36,  0.492,  0.4994,  0.49988,  0.5 

3)  81  elements: 

Ci  =0.15,  0.25,  0.36,  0.43,  0.47,  0.492,  0.4994,  0.49988,  0.5 
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and  additional  mesh  refinement  in  the  corner  if  the  solution  is  singular. 

Also  other  plates  shapes  are  considered  in  the  paper.  We  assume  homogeneous 
isotropic  material  with  Poisson  ratio  v  =  0.3. 

In  Section  2  we  define  the  plate  problem  as  a  three-dimensional  problem 
of  elasticity  and  give  basic  preliminaries. 

In  Section  3  we  define  a  hierarchical  family  of  plate  models. 

Section  4  deals  with  the  problem  of  the  simply  supported  plate.  In  the 
first  part,  4.1,  we  discuss  the  definition  of  the  hard  and  soft  support  for  a 
3-dimensional  formulation.  We  introduce  various  theorems  showing  that  as 
d — >0,  the  difference  between  the  two  supports  disappears  when  measured  in 
the  energy  norm.  A  paradoxical  behavior  of  the  hard  support  shows  that  the 
difference  between  hard  and  soft  support  need  not  necessarily  be  constrained 
to  the  small  neighborhood  of  the  boundary.  It  depends  on  the  structure  of 
the  boundary  of  the  plate.  In  4. 1  we  discuss  the  singular  behavior  of  the 
three-dimensional  solution  in  the  neighborhood  of  the  edges  and  vertices  and 
the  behavior  of  the  moments  and  shear  forces  computed  from  it. 

Section  4.2  presents  detailed  numerical  results  for  the  uniformly  loaded 
square  plate  (in  3  dimensions)  which  serves  in  the  next  section  as  the  exact 
solution.  We  concentrate  on  the  difference  of  moments  ard  shear  forces  for 
the  soft  and  hard  supports  and  descriptions  of  the  boundary  layer. 

Section  4.3  describes  the  family  of  the  dimensionally  reduced  problems 
and  formulates  various  pertinent  theorems.  The  well-known  Reissner-Mindl in 
model  appears  as  one  member  of  this  family.  The  theoretical  behavior  of  the 
solution  of  these  models  in  the  neighborhood  of  the  domain  corners  is 
presented. 

Section  4.4  reports  various  numerical  results.  It  shows  the  accuracy  of 
the  solution  of  various  models  when  compared  with  3-dimensional  solutions  for 
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the  hard  and  soft  support.  The  boundary  layer  and  the  character  of  the 
solution  in  the  neighborhood  of  a  vertex  is  addressed. 

In  Section  4.5,  we  address  the  problem  of  the  L-shaped  domain  for  the 
soft  simple  support.  The  main  emphasis  is  given  to  the  theoretical  and 
numerical  analysis  of  the  solution’s  behavior  in  the  ne’ghborhood  of  the 
vertex  with  the  concave  angle.  The  main  characteristic  difference  in  the 
comparison  with  the  case  of  square  domain  is  shown  .  It  is  shown,  among 
other  things,  that  the  different  models  have  solutions  with  the  different 
strengths  of  the  singularity. 

Section  4.6  analyzes  the  Kirchhoff  model.  It  is  shown  here,  for 
example,  that  the  Kirchhoff  model  essentially  approximates  the  hard  support 
and  the  classical  computation  of  the  reaction  is  aui  attempt  to  get  the  value 
of  the  reaction  for  the  soft  support. 

Section  4.7  summarizes  the  results  presented  in  Section  4. 

Section  5  discusses  some  other  boundary  conditions.  In  Section  5. 1  the 
case  of  square  clamped  plates  is  analyzed.  It  is  shown  aunong  others,  that 
the  Relssner-Mindlin  model  which  captures  well  the  boundary  layer  of  the 
3-dimensional  solution  for  the  simple  support  is  completely  missing  a  strong 
boundary  layer  for  the  clamped  boundary  condition. 

Section  5.2  addresses  then  the  square  plate  which  is  clamped  on  two 
opposite  sides  and  free  on  two  opposite  sides.  Attention  is  given  to  the 
behavior  of  the  solutions  of  various  models  in  the  neighborhood  of  the 
vertices. 

Section  6,  the  final  one,  malces  additional  remarks  and  some  basic 
conclusions. 


4 


2.  Basic  notations  said  preliminaries. 

2 

Let  u  e  IR  be  a  polygonal  domain  with  the  boundary  r  and  Q  =  {x  = 

3  d  d 

(x^X^Xg)  €  R  |  (Xj.Xg)  6  W,  ~2  <  X3  <  n  wil1  k®  called  the  Piate  of> 

thickness  d.  Further  let 

3  d  d 

S  =  {x  €  R  !  (Xj.Xg)  €  r,  <  x2  <  and 

R+  =  {x  e  R3  |  (x1>x2)  e  u,  x3  =  ±|> 

By  the  (exact)  solution  of  the  plate  problem,  we  will  understand  the 
solution  of  the  three-dimensional  linear  elasticity  problem  for  an  isotropic 
homogeneous  material  when  the  equal  normal  load  is  acting  on  R+.  More 
precisely  denoting  T  =  (l^.T^.T^)  the  traction  vector,  we  assume  T  =  T2  = 

l 

0  on  R+  and  T3  =  ^q(x),  x  €  w  on  R+.  On  S  various  homogeneous 
boundary  condition,  which  will  be  specified  later  will  be  considered.  By  u 
=  (u^,  i  =  1,2,3,  and  e  =  {e^},  i,  j  =  1,2,3,  <r  =  {o-^},  i ,  J  =  1,2,3 
we  denote  the  displacement  vector,  the  strain  tensor  and  the  stress  tensor. 

Let  A  =  {a^>,  i,j  =  1,...  ,6  be  the  Hooke’s  law  matrix  (compliance 
matrix)  relating  the  strain  and  stress 


vll' 

‘all 

al2 

al3 

0 

0 

0  ' 

’'ll’ 

°22 

a21 

a22 

a23 

0 

0 

0 

C22 

°33 

= 

a31 

a32 

a33 

0 

0 

0 

C33 

°12 

0 

0 

0 

44 

0 

0 

C12 

°23 

0 

0 

0 

0 

a55 

0 

C23 

/is. 

0 

0 

0 

0 

0 

&66 

C13 

(2.1)  will  be  written  in  the  form 

(2.2)  <r  -  Ac. 
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In  the  case  of  the  isotropic  material  we  have  in  (2.1) 


all  =  a22  =  a33  =  X  +  =  (1  i  v)(l}-  2v) 

(2'3)  al2  =  a13  =  a23  =  *21  ~  a31  =  a32  ~  X  =  (  l"  +  u )  ( 1  -  2u ) 
a44  =  a55  =  a66  =  2M  =  ITU  * 


A,p  are  the  Lam6  constants,  E  is  the  modulus  of  elasticity,  v  is  the 
Poisson’s  ratio,  and 


A  = 


vE 

(1  +  v)(l  -  2v) 


E 

M  2(1 + v) 


Further  we  let 


(2.4) 


A,  , 
e  (u) 


f  (I  dx 

J  i  1=1 
n 


r  lleil  +  <r22®22  +  <r33e33  +  2(rl2el2  +  2°23e23  +  2<rl3el3 


dx 


A 

be  the  strain  energy  expression.  We  use  the  notation  e  (U)  to  emphasize 
that  (2.4)  is  based  on  A  given  by  (2.1).  By  the  total  energy  we  denote 
as  usual 


(2.5) 


GA(u)  =  e^(u)  -  Q(u) 


where 


(2.6) 


Q(u)  =  |  2 (U3 ( x i ’ x2 ’ 2 ^  +  U3(xi’x2,-2^)dx 


u> 


The  exact  solution  u  of  the  plate  problem  is  the  minimizer  of  the  total 
energy  over  the  set  of  functions  Jf(fl)  c  (H1^))3  where  Jf(fl)  constrains 
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1  3 

(H  (fl))  on  S  (not  on  R+).  The  boundary  conditions  of  the  plate  problem 
are  then  uniquely  characterized  by  the  set  H(Q). 


3.  The  hierarchic  family  of  plate  models. 

By  the  plate  model  we  mean  a  two-dimensional  boundary  value  problem 
which  approximates  the  solution  of  the  three-dimensional  plate  problem 
(defined  in  Section  2).  By  the  hierarchic  family  of  plate  models  we 
understand  a  sequence  of  models  which  solutions  converge  to  the  exact 
(three-dimensional)  solution  of  the  plate  model  and  any  model  of  the  sequence 
converges  (after  scaling)  as  d — >0  to  the  same  limit. 

The  plate  model  from  the  hierarchic  family  is  defined  as  the  minimizer 
u(n)  of  the  quadratic  functional 

(3.1)  GB(u)  *  eB(u)  -  Q(u) 

over  the  set  K(n)  c  H(fl)  of  the  functions  of  the  form 


(3.2) 


(n) ,  . 
u.  (x) 


ni 


=  1 
J=0 


u[|J)(x1,x2)(x3)  . 


i  =  1,2,3 


n  =  (n.  ,n0,n„)  and  B  is  certain  compliance  matrix  depending  on  n,  which 

i  C  O 

could  be  different  from  A.  The  solution  based  on  (3.1)  and  (3.2)  will  be 
called  the  solution  of  the  n-model. 

Because  of  the  assumptions  about  the  symmetry  of  the  load  we  made  in 
Section  2  we  can  assume 


UU  =  U2J  ■  ° 


for  j  even 


and 


“3J  =  ° 


for  J  odd. 
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Often,  instead  of  (3.2), 


n‘ 

(3.2a)  u[n,(x)  =  Y.  “u'VVjMj 

J=0 


is  used.  Here  is  the  Legendre  polynomial  of  degree  J.  The  solution 

u(n)  using  (3.2)  or  (3.2a)  is  the  same  because  the  span  is  the  same.  Form 
(3.2a)  is  preferable  especially  when  adaptive  approaches  are  used. 

Finally  by  c  Jt  we  denote  the  set  of  functions  of  the  form 

ui(xl'x2’x3)  =  “  diii  X3*  1  =  1,2 

(3.3) 

U3^X1’X2’X3^  =  U3^X1,X2^ 


with  u3(x1,x2)  €  H^u).  Function  uk  is  then  the  minimizer  of  the 
quadratic  functional  (3.1)  over  3^  (when  using  matrix  B  which  will  be 
specified  later). 

It  is  obvious  that  the  set  J((n)  is  dense  in  H(n)  as  n 


‘i 


l  = 


1,2,3  and  hence  ^u^n^ 


>u  as  n. 


>«,  i  =  1,2,3  with  the  convergence 


in  the  energy  norms.  By  index  A  in  *u^  we  indicate  that  the  matrix  A 
in  (3.1)  is  used. 

There  is  a  vast  amount  of  literature  devoted  to  the  derivation  of 
various  plate  models.  We  refer  here  to  [1,3,4]  and  references  therein. 


4.  The  problem  of  the  simply  supported,  uniformly  loaded  (q(x)=l)  plate. 

In  this  section  we  still  analyze  various  features  of  the  plate  problem 
when  the  simple  support,  which  is  the  typical  one  in  engineering  analysis,  is 
Imposed  at  the  boundary.  In  Section  5  we  will  address  the  plate  problem  with 
some  other  boundary  conditions  and  discuss  the  properties  of  its  solution  in 
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comparison  with  the  solution  of  simply  supported  plates. 

4.1  The  (three-dimensional )  plate  problem  and  the  basic  properties  of  its 
solution. 

In  this  section  we  will  address  the  three-dimensional  plate  problem  and 
the  properties  of  ltr  solution.  The  solution  will  serve  as  the  basis  for  the 
assessment  of  the  accuracy  of  various  plate  models. 

The  simple  support  is  an  idealization  which  has  no  standard  definition 
in  the  three-dimensional  setting.  We  will  consider  here  only  two  different 
types  of  simple  s-'oports: 

a)  the  soft  simple  support 
0)  the  hard  simple  support. 

The  soft  simple  support  is  characterized  by  the  set  3^(0)  =  {u  e 
1  3 

CH  (n))  |  u3  =  0  on  S>  and  the  hard  simple  support  is  characterized  by  the 

set  K^(n)  =  {u  e  ( H1  ( n) ) 3 1  u3  =  0,  u^  =  0  on  S>.  Here  by  ufc  we  denote 

the  displacement  in  the  tangential  direction  cn  T,  i.e.  u^  =  u^tj  +  u^^ 

where  t  =  (t^,t2)  is  the  unit  tangent  to  T.  In  the  vertices  of  w  (u  is 

assumed  to  be  a  polynomial  domain)  no  constraint  condition  is  prescribed. 

1  3 

Because  u  e  (H  C*D)  the  constraint  is  interpreted  in  the  usual  trace 
sense. 

The  basic  physical  interpretation  of  both  boundary  conditions  is 
obvious.  In  the  case  of  the  soft  support,  no  tangential  shear  stresses  aire 
present,  while  in  the  case  of  hard  support  the  tangential  stresses  lead  to 
the  twist  moment  reaction. 

Let  us  first  discuss  the  differences  between  the  twe  mentioned  cases  of 
simple  support  for  d  small. 

S  H 

Theorem  4.1.  Let  u,  and  u"  be  the  solution  of  the  plate  problem  with  the 

a  n 

thickness  d  for  the  soft  and  hard  support  cases,  respectively.  Then 
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(4.  1) 


A,  S  H, 
e  (ud  -  ud) 

A,  S, 
e  ( ud ) 


ELS  d - >0 


The  proof  follows  by  arguments  used  in  [2]  [6]. 


□ 


1/2 

Because  If5  d  we  have  eA(u^)  £  eA(u|j)  and  [/(u^-u^)]  = 

,  ,1/2  d  d  l  d  d  J 

e  (ud)-e  (u^)  •  Theorem  4.1  shows  that  the  difference  between  both 

supports  measured  in  the  (relative)  energy  norm  converges  to  zero  as  d - >0. 

It  is  necessary  to  underline  that  this  holds  for  the  energy  norm  and  not  for 
all  other  norms  els  will  be  seen  later.  The  main  difference  between  the  two 
supports  is  the  behavior  of  the  solution  in  the  neighborhood  of  S.  The  size 
of  this  neighborhood  depends  on  the  relation  between  the  plate  thickness  d 
and  the  smoothness  of  T.  The  following  example  shows  that. 

Let  us  consider  to  be  a  regular  m-gon  inscribed  in  the  unit 

circle  C  as  shown  in  Figure  4.1. 


Figure  4.1.  Scheme  of  the  m-gon  u  ,  (m  =  6). 

m 
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H  S 

Denote  by  u,  and  u.  ,  respectively,  the  (three-dimensional)  solution 
a,  m  a, m 

H  S 

for  the  hard  and  soft  support  for  the  m-eck  w  .  By  u ,  and  u ,  we 

m  d,c  d,c 


denote  the  solution  on  C  x 


[r  2] 


where  on  T  x 
c 


('2-  2] 


the  hard  and 


soft  supports,  respectively,  are  imposed.  The  solution  is  defined 
analogously  as  the  minimizer  of  the  energy  over  ^  or  We  then  have 


Theorem  4.2.  For  any  0  <  d  <  d^, 


lim 
m — x» 


S  S 

u  ,  — » u  , 

d,m  d,c 


lim 
m — x» 


1,  ^  U, 

d,c  d,c 


The  convergence  is  understood  in  the  energy  norm  on  12  .  For  the  proof 


see 

[6].  □ 


It  is  possible  to  compute  the  value 

^(i/)  =  lim  lim  |u».d(0’°>0)  ~  u5,d.L0,0,0),l  t 
d — jO  m— ( uB  t  d  ( 0 , 0 , 0 )  | 

£(0.3)  =  0.264 

(We  have  lim  |uS  ,(0,0,0) |  >  lim  |u^  ,(0,0,0|) 
m,  d  m,d 


Theorem  4.2  shows  that  when  m  is  large  with  respect  to  d  *  then  the 
neighborhood  of  S  where  there  is  significant  difference  between  hard  and 
soft  support  solution  can  "fill"  the  entire  fl. 

The  boundary  of  the  plate  12  has  edges  and  vertices.  Local  behavior  of 
the  solution  of  the  plate  problem  in  the  neighborhood  of  these  edges  and 
vertices  is  well-known.  See  for  example  [7]  and  [8]. 
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Let  us  consider  .'irst  the  horizontal  singularity.  To  this  end  let  us 
assume  that  1^  is  a  side  of  the  polygon  w  and  =  (x^Xg  |  I x1 1  <  a, 
x2  =  0>.  Then  the  (horizontal)  edge  is  Ej  =  (Xj.Xg.Xgl  (x^Xg)  e  X3  = 
^}.  In  the  neighborhood  of  E^  we  have 

UjCXj.Xg.Xg)  =  +  smoother  term 

(4.2)  u2(x1>x2,x3)  =  C(x1)r^2(0)  +  smoother  term 

u3(xlfx2,x3)  =  C(x1)r^03(0)  +  smoother  term 


where  (r^,0)  denotes  the  polar  coordinates  in  the  plane  x2,x3  with  the 
origin  in  E^.  Functions  ^2  and  03  are  analytic  in  0  and  function 
C(x^)  is  smooth  on  (open)  interval  -a  <  x^  <  a.  In  (4.2),  the  smoother 
term  „ieans  a  function  which  is  smoother  in  the  neighborhood  of  Ej  than  the 
functions  r^^fo),  1  =2,3.  For  the  hard  and  soft  support  we  have  X  = 
0.71117  (independently  of  the  Poisson’s  ratio  v).  Because  X  <  1,  the 
stresses  are  unbounded  in  the  neighborhood  of  E^.  Nevertheless  the  moments 


(4.3) 


.d/2 

VvV  "J  ■rijx3dv 

-d/2 


i.J,  =  1,2 


and  the  shear  forces 


(4.4) 


Q3,i(Xl 


d/2 

x2)  *  J  °'ljdX3' 

-d/2 


i  =  1,2 


are  bounded  in  the  neighborhood  of  (except  possibly  the  vertices  of 

V- 

Consider  now  the  vertical  edge  E^.  E^  =  (x^  *  0,  x2  =  0,  | X3 |  < 
when  we  assumed  that  the  vertex  of  u  lies  in  the  origin.  We  will  assume 
that  the  Internal  angle  at  this  vertex  is  a,  0  <  a  £  2n.  In  the 


12 


neighborhood  of  Eg  we  have 

ul^xrx2,x3^  =  Cl^x3^r21<^l^9^  +  smoother  term 
(4.5)  UgfXj.Xg.Xg)  =  CjCXgJr^^gCe)  +  smoother  term 

Ugfx^Xg.Xg)  =  C2(x3)rg2^(0)  +  smoother  term 

Here  (rg,0)  31,6  the  P°lar  coordinates  in  the  x^,Xg  Plane.  ^(0)  and 

and  Cg(x 

with  possible  singular  behavior  at  ±^. 

Coefficients  A^  and  Ag  and  functions  0  depend  on  a  and  the 

type  of  boundary  conditions.  Coefficient  A^  can  be  complex.  Then  there  is 
a  pair  of  conjugate  coefficients  because  (4.5)  has  to  be  real. 

Ul(xl,x2’x3^  “  Cl^x3^r2C^1  sin(ymIm A1lgr)  0X(0) 

+  c1(x3^r2e^1  cos(ymIm  AJlgr)  ^1(0) 

+  smoother  terms 


^(0)  are  analytic  in  0,  and  C^(Xg) 


)  are  smooth  on 


[■*•  2] 


and  analogously  UgfXj, Xg, Xg) . 

The  coefficients  A^,Ag  can  ^  inte8ers  or  they  can  be  of  multiplicity 

>1  for  some  angles  a.  Then  these  angles  are  exceptional  and  in  general  the 
A  AS 

term  r  has  to  be  replaced  by  r  (lgr)  ,  S  an  integer.  Only  a  finite 
number  of  exceptional  angles  exist.  In  the  neighborhood  of  these  exceptional 
angles  the  singular  behavior  is  not  continuous  with  respect  to  the  a.  We 
will  not  address  here  these  cases  although  they  will  also  occur  in  some  of 
our  examples. 

In  Table  4.1  we  show  the  values  of  A^  and  Ag  in  (4.5)  for  hard  and 
soft  supports  for  various  angles  a.  They  are  independent  of  the  Poisson’s 
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ratio  v.  We  see  that  the  angles  a  =  30°  45°  90°  are  exceptional  because 
Ag  is  sm  integer.  The  angle  360°  leads  to  the  multiplicity  two.  For  a  = 
30-135  the  coefficient  A^  in  the  table  for  the  soft  support  case  are 
complex.  For  more  see  [9]. 


Table  4.1  The  coefficients  and  A2  in  (4.5)  for  the  hard  and  soft 

support . 


a 

*1 

a2 

HARD 

SOFT 

HARD 

SOFT 

30 

3.4846 

8.0630 

+ 

i  4.2028 

6.0000 

6.0000 

45 

2.4129 

5.3905 

+ 

i  2.7204 

4. 0000 

4.0000 

90 

1.4208 

2. 7396 
+ 

i  1.1190 

2.0000 

2.0000 

120 

1.2048 

2.0941 

+ 

i  0.6046 

1.5000 

1.5000 

135 

1. 1368 

1.8853 

+ 

i  0.3606 

1.3333 

1.3333 

150 

1.0832 

1.5339 

1.2000 

1.2000 

225 

0.7263 

0.6736 

0.8000 

0.8000 

270 

0.5951 

0.5445 

0.6667 

0.6667 

315 

0.5330 

0.5050 

0.5714 

0.5714 

360 

0.5000 

0.5000 

0.5000 

0.5000 

We  see  that  for  some  angles  we  have  A  <  1  and  hence  in  these  cases  the 
stresses  are  unbounded  as  — »0.  We  have  then 
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(4.6) 

M1J  ■  Clr21'‘*lJle> 

+  smoother  terms, 

i.J  =  1,2 

(4.7) 

=  C2r2rWe> 

+  smoother  terms, 

J  =  1.2 

We  mention  that  the  expressions  (4.6)  and  (4.7)  do  not  follow  directly  from 
(4.5)  because  the  solution  has  other  singularities  in  the  vertex  and  also 
C1(x3)  in  (4.5)  can  be  singular  for  x3  =  .  Nevertheless  these 

singularities  are  weaker  and  do  not  influence  the  form  (4.6)  (4.7). 

The  expressions  (4.6)  and  (4.7)  show  that  the  moments  and  shear  forces 
in  the  neighborhood  of  the  vertex  of  u  can  be  unbounded.  For  example,  this 

O 

happens  when  u  is  an  L-shaped  domain  and  a  *  270  .  The  shear  forces  have 
the  same  strength  of  the  singularity  for  the  hard  and  soft  support,  while  the 
singularity  of  the  moments  is  different  for  the  hard  and  soft  support.  Fur¬ 
thermore,  the  singular  behavior  is  different  for  the  moments  and  shear 
forces.  The  mentioned  behavior  described  by  (4.6)  and  (4.7)  is  valid  only 
for  r  «  d  as  will  be  seen  later  in  Section  4.5. 

It  is  interesting  to  mention  that  the  differences  in  the  singular 
behavior  of  the  hard  and  soft  supports  do  not  explain  Theorem  4.2,  whose 
proof  is  based  on  completely  different  principles. 

4.2  The  (three-dimensional)  problem  of  a  simply  supported^  uniformly  loaded 


square  plate.  Numerical  results. 


In  this  section  we  will  present  numerical  results  of  the  analysis  of  the 
square  plate. 

Let  u  =  {Xj.Xg  1  <  °-5>  I X2 1  <  °-5^  and  let  d  =  f(5o  511101 

consider  the  three-dimensional  problem  with  v  -  0.3.  In  this  case  a  =  90° 
and  hence  from  Table  4. 1  we  conclude  that  the  moments  and  shear  forces  are 
bounded  and  are  sufficiently  smooth  up  to  the  boundary. 


15 


The  main  question  we  will  consider  in  this  section  is  the  difference 


between  the  solution  when  the  hard  and  soft  simple  support  is  prescribed 
on  S.  The  solution  is  obviously  symmetric  with  respect  to  the  axes  x^^  and 
Xg  and  hence  oriy  a  quarter  of  the  plate  will  be  considered. 

First  we  consider  the  behavior  of  the  twist  moment  for 

=0.4,  Xj  =  0.5  as  a  function  of  the  variable  Xg.  Denoting  M^2  and 
respectively,  the  twist  moments  for  the  hard  and  soft  support,  we 
have  =  0  while  This  indicates  that  the 

boundary  layer  has  to  be  present  because  of  Theorem  4. 1  which  shows  that  the 
difference  of  both  solutions  converges  to  zero  as  d — >0. 

In  Table  4.2  we  show  the  values  of  the  twist  moment  Mi2^Xl,X2^‘  We 
clearly  see  that  the  boundary  layer  is  of  order  d. 

Table  4.2  The  twist  moment  M^CXj.Xg)  x^  =  0.4,  0.5  for  the  hard  and 
soft  support. 


X2 

X1  S 

=  0.4 

X1  ' 

=  0.5 

HARO 

SOFT 

HARD 

SOFT 

0.0 

0. 

0. 

0. 

0 

0.02368 

0.0018 

0.0019 

0.0020 

0 

0. 11842 

0.0092 

0. 0092 

0.0097 

0 

0.21316 

0.0162 

0.0163 

0.0172 

0 

0.45000 

0.0289 

0.0292 

0.0315 

0 

0.48079 

0.0294 

0.0297 

0.0323 

0 

0.49026 

0.0295 

0.0285 

0.0324 

0 

0.49500 

0.0295 

0.0237 

0.0325 

0 

0.49713 

0.0295 

0.0179 

0.0325 

0 

0.49903 

0.0295 

0.0080 

0.0325 

0 

0.49950 

0.0295 

0.0045 

0.0325 

0 

0. 50000 

0. 0295 

0. 

0.0325 

0 

Remark.  The  exact  solution  is  not  known.  Nevertheless  the  data  we  report 
are  reliable.  They  were  computed  by  refined  meshes,  and  high  degree  elements 
with  am  amalysis  of  the  accuracy. 

Figure  4.2  shows  M^2(0. 4,  Xg)  with  maxM^2(0.4,  Xg)  =  0.0295  and 
M^2(0.4,  x2)  with  maxM^2(0.4,  x2)  =  0.0297. 


<N 


Figure  4.2.  The  moments  M^2(0.4,  x2)  and  M^2((0.4,  x2)  for  the 

3-dimensional  problem. 

The  shear  forces  Q^tx^Xg),  J  ”  1,2  show  the  similar  boundary  layer  we 
mention  that  0^2^xr  “  0  while  (^(Xj,  *  °- 
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Table  4.3.  The  shear  force  Q^jCx^.Xg),  J  =  1,2,  x,  =  0.4,  0.5  for  the 


hard  and  soft  support. 


HARD 

SOFT 

HARD 

SOFT 

HARD 

SOFT 

0.0 

0.2460 

0.2446 

0. 

0.0 

0.3373 

0.4206 

0.02368 

0.2451 

0.2445 

0.0041 

0.0041 

0.3373 

0.4206 

0. 11842 

0.2347 

0.2346 

0.0215 

0.0214 

0.3260 

0.4080 

0.21316 

0.2086 

0.2086 

0.0386 

0.0386 

0.2973 

0.3740 

0.45000 

0.0526 

0.0503 

0. 1274 

0. 1274 

0. 1073 

0. 1326 

0.48079 

0.0207 

-0.0046 

0. 1526 

0. 1526 

0.0534 

-0. 1334 

0.49026 

0.0107 

-0.4193 

0. 1606 

0. 1633 

0.0306 

-0.8126 

0.49500 

0.0053 

-1.9206 

0. 1660 

0. 1760 

0.0184 

-2. 1746 

0.49713 

0.0030 

-3.7753 

0. 1680 

0. 1880 

0.0133 

-3. 5240 

0.49903 

0.0012 

-6.9740 

0. 1701 

0.2066 

0.0038 

-5.4623 

0. 49950 

0. 0005 

-8. 1673 

0. 1706 

0.2133 

0.0019 

-6.0693 

0.50000 

0.0 

-9.7186 

0. 1713 

0.2213 

0.0 

-6.7933 

0 32 


SOFT 


0.0 

0.0  -0.6433 

0.0  -3.1986 

0.0  -5.6633 

0.0  -10.4513 

-10.6353 
-9.9046 
-8.6002 
-7.7056 
-6.9546 
-6.8382 
-6.7933 


In  Figures  4.3  and  4.4  we  show  Q^(0.5,  x,,)  (=  reaction)  and 


023(0.5,  Xg)  for  the  soft  support. 
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Figure  4.3.  The  shear  force  0^(0.  5,  x^)  for  the  soft  support 
(three-dimensional  problem). 


Figure  4.4.  The  shear  force  (^(O.S,  x^) 


for  the  soft  support 


(three-dimensional  problem). 
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Tables  4.2,  4.3  and  Figures  4.2,  4.3,  and  4.4  show  that  the  difference 
between  hard  and  soft  support  resides  in  the  boundary  layer  behavior.  The 
hard  support  does  not  have  a  boundary  layer  and  can  serve  as  a  smooth  exten¬ 
sion  to  the  boundary.  To  understand  the  boundary  layer  more  precisely  we  let 


M?2(xi,x2)  ~Mf2(x1,x2) 
M?2(xlt  0.5) 


exp 


(xl*x2) 


(0.5  -  x2)l 

J 


and 


Q3l(Xl,X2)  -  q!i(xi,x2) 

Q3,i(xi,  0.5) 


exp  [-P31(xvx2^°'--d~  Xz)) 


Table  4.4  shows  some  values  of  &i2^X\,X2^  ^31^X1,X2^’ 


Table  4.4  The  functions  P^CXj.Xg)  and  P31(xlfx2). 


X2 

01 2 ( Xi , X2 ) 

@3 1 ( Xj , X2 ) 

Xi  =  0.4 

X!  *  0.49 

xa  =  0.4 

X!  =  0.49 

X!  =  0.5 

0.49263 

3.34 

2.74 

3.22 

3.72 

2.  17 

0.49500 

3.28 

3.  10 

3.24 

3.53 

2.27 

0.49666 

3.26 

3.  18 

3.28 

3.53 

2.28 

0.49878 

3.27 

3.30 

3.40 

3.62 

2.23 

0.49950 

3.31 

3.34 

3.51 

3.73 

2.22 

0.49989 

3.46 

3.49 

3.62 

3.84 

2.22 

0.49991 

3.52 

3.54 

3.63 

3.84 

2.21 

0. 49994 

3.62 

3.64 

3.63 

3.84 

2.21 

In  Figure  4.5  we  show  the  function  0^  and  in  Figure  4.6  we  show  the  0^. 
We  see  that  the  03^  and  032  increases  when  x2~ >0.5.  This  is  typical 
for  the  three-dimensional  formulation  where  essentially  infinite  number  of 
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boundary  layers  are  present  but  here  only  the  first  one  is  visible.  Also 
typical  is  the  discontinuity  of  031 (x^ , Xg)  when  x^^ — >0.5.  We  see  that 
for  x1  =  0.5  the  value  of  031  is  significantly  smaller  than  for  x1  =0.49. 


0.490  0.492  0.494  0.496  0.498  0.500 


Figure  4.5.  The  function  ^(Xj.Xg)  for  Xj  =  0.4,  0.49,  0.5 
(three-dimensional  formulation). 


0.490  0.492  0.494  0.496  0.498  0.500 


Figure  4.6.  The  function  ^j2^X1'X2^  ^or  X1  =  ^ 

(three-dimensional  formulation). 
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The  shear  force  0^(0. 5,  x^)  is  the  reaction  of  the  plate.  We  see  from 
Fig.  4.3  and  Table  4.3  that  for  the  soft  support  this  reaction  is  negative  in 
the  neighborhood  of  the  vertex  x2  *  0.5.  This  relates  to  a  well-known  effect 
in  the  Kirchhoff  theory  (see  Section  4.6).  It  is  Interesting  to  define 

3?  «  J  Qgj(0. 5,  x2)  dx2 

I 


where  I  =  (0  <  x2  <  0.5  |  0^(0.  5,  Xg)  i  0).  We  get  in  our  case  3?  =  0.0290. 

The  data  we  present  here  are  for  d  =  •  Q3  ^(O.S,  0.5)d — » 

^0.5,  0.5),  as  d — >0  where  ^  is  a  certain  value  which  we  expect 
to  be  approximately  0.068,  although  the  authors  do  not  know  whether  a 
rigorous  proof  of  it  exists.  Analogously  X — as  d — »0.  We  expect 
that  W  *  0.029  which  is  different  from  the  value  0.0325  arising  from  the 
Kirchhoff  model  (see  Section  4.6). 

The  boundary  layer  is  also  visible  in  the  energy  of  the  solution.  We 
have  in  our  case  d  = 


'  A, 
e  (u  ) 


A.  S. 
c  (u  ) 


1/2 


=  0.895 


and  hence  the  difference  between  hard  and  soft  support  measured  in  the  energy 
norm  is  954.  Of  course  as  d — >0  the  energy  norm  of  this  difference 
converges  to  zero  due  to  Theorem  4.1.  We  discussed  only  the  square  plate 
here.  If  the  plate  is  rhombic  then  the  difference  between  the  hard  and  soft 
support  Increases. 


4.3  The  dimensionally  reduced  models. 


We  will  consider  here  the  n-model  Introduced  in  Section  3  which  is  based 

g 

on  the  minimization  of  the  quadratic  functional  G  (u)  given  in  (3.1)  over 
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the  set  of  all  functions  of  the  form  (3.2)  constrained  in  the  space  J((Q). 

We  have  now: 

Theorem  4.3.  Let  n  =  (n^.^.n^),  n^  i  1,  n^  i  1,  n^  2  2.  Denoting  by 
V,n)  the  plate  model  based  on  (3.1)  (3.2)  with  B  =  A  (A  given  by  (2.1) 
(2.3))  then  for  both  types  of  simple  support  we  have 

A, A  n  , 

(4.8)  c  (  ud  -  Uj2  - >  o  as  d — >0 

cA(ud) 

The  proof  follows  by  using  arguments  in  [6].  □ 


The  assumption  that  n3  2  2  in  Theorem  4.3  is  essential.  If  n^  =  0 
and  v  =  0  then  (4.8)  still  holds  but  when  v  >  0  (4.8)  does  not  hold. 
This  means  that  the  case  n^  =  0  with  the  matrix  A  cannot  be  used  except 
when  v  =  0.  Hence  the  following  problem  arises:  Find  a  class  2  of 
matrices  B  which  when  used  instead  of  A  leads  to 


(4.8a) 

Here  by  Bu^ 

B  inste  ad  of 
The  class 


A,B  n  , 
e  (  ud  -  Ud) 


cA(u<i) 


1/2 


->  0 


as 


we  denoted  the  minimizer  of  (3.1)  with 


d — > 0 

n  =  (1, 1,0) 


and  matrix 


A. 

2  which  guarantees  (4.8a)  is  the  class  of  matrices  (4.9). 
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(4.  9)  B  = 


l-i> 

E 

l-i/ 

0 


uE 

1-/ 

E 

l-i/£ 

0 


0 

E 


1-y 

0 


E 

1+v 


0 

EH 

1+1/ 

0 


EH 

1+w 


Here  E  and  v  are  the  modulus  of  elasticity  and  the  Poisson’s  ratio  and 
H  >  0  is  an  arbitrary  parameter.  We  have. 

Theorem  4.4.  Let  n  *  (1,1,0)  and  let  Bun  be  the  minimizer  Of  (3.1)  (3.2) 
with  the  matrix  B  given  by  (4.9).  Then  for  both  types  of  simple  support  we 
have 


(4.9) 


f  A,  B  n,1 
g  (ud  -  ud) 

.  eA(ud) 


1/2 

- >  0 


as  d  — >0 


The  proof  follows  by  the  arguments  used  in  [6]. 


□ 


The  coefficient  H  in  B  does  not  influence  the  validity  of  (4.9)  but 

it  influences  the  size  of  the  difference  BUj-u..  The  solution  Bu^  (n  * 

d  d  d 

(1,1,0))  is  the  famous  Reissner-Mindlin  solution.  The  recommendation  for 
the  selection  of  H  was  made  by  various  authors.  In  [10]  the  value  5/6  is 
recommended.  In  the  connection  with  wave  propagation  in  [11],  suggested 
values  range  almost  from  0.76  for  v  =  0  to  0.91  for  v  =  0.5.  Another 
recommendation  is  made  e.g.  in  [12].  Jf(0)  =  0.833,  X(0.33)  =  0.850,  H(0.5) 

*  0.87.  We  will  address  this  question  later  in  the  numerical  way. 

Let  us  mention  that  for  any  n-model  with  n^  £  2,  using  the  matrix  B 
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will  not  lead  to  (4.9). 


The  n-model  reproduces  Theorem  4.2.  We  have, 

Theorem  4.5.  Let  £  2  with  using  A  or  n3  =  0  using  B. 


Then  for  any 


dSd, 


11m  u 


m 


(n)S 

ld,n 


u 


(n)S 
d,  c 


Urn  u'n)H  - ,  S'n)  «  u'n) 

d,n  d,c  d,c 


m 


For  the  proof,  again  see  (6]. 


□ 


Using  matrix  B  and  (3.3)  together  with  (3.1)  (3.2)  we  obtain  the 

Kirchhoff  solution  u,  and 

a 

Theorem  4.6.  We  have 


(4.10) 


A,  k, 

e  (ud  -  ud) 


eA(ud) 


1/2 


0 


as  d 


See  [6]  for  proof. 


□ 


The  Kirchhoff  model  cannot  distinguish  between  the  soft  and  hard 
support  (see  Section  4.6  for  more)  and  the  second  part  of  Theorem  4.2  holds 
for  the  Kirchhoff  model  too.  Various  paradoxical  properties  of  Kirchhoff 
model  for  simple  support  were  discussed  in  [6,13,14]. 

The  n-model  solution  has  a  boundary  layer.  In  [15]  and  [16]  the  boundary 
layer  of  the  Reissner-Mindlin  solution  was  rigorously  analyzed  in  detail  when 
u  has  smooth  boundary.  The  rigorous  analysis  of  the  boundary  layer  structure 
of  the  solution  for  a  domain  u  with  piecewise  smooth  boundary  is  not  avail¬ 
able.  For  additional  analyses  see  also  [17].  The  solution  from  the  Kirchhoff 
model  has  no  boundary  layer. 

From  (4.9)  and  (4.10)  it  follows  that  the  Reissner-Mindlin  solution  con- 
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verges  in  energy  to  the  Kirchhoff  model.  Nevertheless  we  underline  that  this 
is  only  in  the  energy.  At  the  boundary  the  convergence  does  not  occur.  Using 
the  analytical  results  in  [15]  for  infinite  half-plane  this  can  be  easily 
analytically  demonstrated. 

In  Section  4.2  we  analyzed  the  singular  behavior  of  the  plate  solution  in 
the  neighborhood  of  the  corners  of  the  domain.  We  especially  addressed  in 
(4.6)  (4.7)  the  singularity  of  and  as  functions  of  internal  angle 

a.  The  Relssner-Mindlin  model,  n-model,  and  Kirchhoff  model  solution  also 
lead  to  the  singular  behavior  of  and  Q^j.  The  forms  (4.6)  (4.7)  are 

also  valid  here,  i.e.  we  have 


C  r  1  1 

Ulr2 


%  (n) 

C  r  2 
L2r2 


(0)  +  smoother  terms. 


^(0)  +  smoother  terms. 


i.j  =  1.2 

J  -  1,2 


In  general  A^  depends  on  n,  i.e.  on  the  model,  and  they  are  different 
for  different  models  and  for  different  types  of  boundary  conditions.  In  Table 
4.1  we  gave  the  values  of  Aj,  i  *  1,2  for  the  three-dimensional  problem. 

In  Table  4.5  we  list  the  coefficients  A^  for  the  n-models  and  Kirchhoff 
model  for  the  soft  support  and  hard  support  analogously  to  Table  4.1.  We 
mention  that  the  Kirchhoff  model  does  not  distinguish  between  the  hand  and 
soft  support. 


Table  4.5  The  coefficients  Aj  and  Ag  for  the  hand  and  soft  for  vanious 
models. 
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*1 

X2 

a 

HARD 

SOFT 

HARD  SOFT 

RM 

(1,1,2) 

3  dim 

RM,  (1,  1,2) 
3  dim 

K 

RM,  (1,  1,2) 
3  dim 

K 

30 

3.7972 

3.4846 

8.0630 

+ 

i  4.2028 

5.0000 

6.0000 

4.0000 

45 

2.6083 

2.4129 

5. 3905 
+ 

i  2.7204 

3.0000 

4. 0000 

2.0000 

90 

1.4905 

1.4208 

2.7396 

+ 

i  1.1190 

1.0000 

2.0000 

0.0000 

120 

1.2404 

1.2048 

2.0941 

+ 

1  0.6046 

0.5000 

1 . 5000 

-0. 5000 

135 

1. 1609 

1. 1368 

1.8853 

+ 

i  0.3606 

0.3333 

1.3333 

-0.6666 

150 

1.0978 

1.0832 

1.5339 

0.2000 

1 . 2000 

-0.8000 

225 

0.7354 

0.7263 

0.6736 

0.2000 

0.8000 

-0.8000 

270 

0.6040 

0.5951 

0.5445 

0.3333 

0.6667 

-0.6667 

315 

0. 5379 

0. 5330 

0. 5050 

0. 1428 

0.5741 

-0.8572 

360-0 

0.5000 

0.50000 

0.5000 

0.0000 

0. 5000 

-1.0000 

For  a  =  360°  we  list  the  limiting  value  a  =  360-0  which  is  not  equal 
to  the  value  a  =  360°.  We  see  that  for  the  hard  support  the  Reissner-Mindlin 
model  and  (1,1,2)  model  give  slightly  different  values  while  the  Kirchhoff 
model  yields  very  different  singularity  coefficients.  For  the  formula  for 
these  values  we  refer  to  [9]. 

4.4.  The  n-model  of  simply  supported  uniformly  loaded  square  plate. 

Numerical  results. 

In  this  section  we  consider  the  same  problem  as  in  Section  4.2  namely 
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the  square  plate  d  =  ,  v  =  0.3  and  soft  and  hard  simply  support. 

7 

Further  we  will  assume  that  E  =  10  . 

As  we  mention  in  Section  4.3  the  Reissner-Mindlin  model  (n  =  (1,1,0)) 
depends  on  the  value  of  the  parameter  Jf.  Let  us  first  show  the  dependence 
of  the  energy  e(Jf)  (of  the  quarter  of  the  plate)  on  the  value  of  Jf  for 
various  models.  Denote  further  the  quantity 

1/2 

The  values  of  e(J 0  and  e  for  the  soft  support  is  shown  in  Table  4.6. 

We  see  that  by  using  different  values  of  the  parameter  Jf,  the  energy  of 
the  RM  model  can  be  over  or  under  the  true  3-dimensional  energy.  The  energy 
in  the  case  (1,1,2)  is  smaller  than  the  energy  of  the  RM  solution  for  Jf  =  1. 
If  v  -  0  then  for  Jf  =  1  we  get  A  *  B  and  the  energy  of  (1,1,2)  model 
would  be  larger  than  that  of  the  RM  model.  It  seems  remarkable  that  the  per- 

5 

formance  of  RM  model  (for  this  particular  case)  is  better  for  all  S  Jf  <  1 

b 

than  that  of  the  model  (1,1,2)  and  Jf  =  1.  This  conclusion  is  not  valid  in 
general  for  all  boundary  conditions.  The  model  (3,3,4)  (Jf  *  1)  gives  then 
the  better  results  than  the  RM  model. 


|e(Jf)  -  e(3  dim)  I 
I e(3  dim)  | 
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Table  4.6  Energy  e(JO  as  a  function  of  H  and  different  plate  models  for 


Table 
support . 

Table  4.7 


the  soft  support 


H 

cm  x  103 

e% 

Model 

1 

0.  234563 

2.69 

RM 

0.91 

0.  234674 

2.26 

RM 

0.87 

0.234729 

0.41 

RM 

0.845333 

0.  234765 

1.  16 

RM 

5/6 

0. 234783 

1.46 

RM 

1 

0 . 234528 

2.95 

(1,1,2) 

1 

0. 234731 

0.29 

(3,3,4) 

1 

0 . 234732 

0.01 

(5,5,6) 

0.87 

0. 234693 

1.3  ' 

(1,1,2) 

0.87 

0.234913 

2.76 

(3,3,4) 

0.87 

0.234914 

2.78 

(5,5,6) 

Klrchhoff 

0.232392 

9.98 

K 

3  dim 

0.234735 

— 

— 

4.6  has  addressed  the  soft  support.  Table  4.7  addresses  the  hard 


The  energy  for  the  various  models  for  the  hard  support. 


H 

e  x  103 

cY. 

Model 

0.87 

0.232524 

1.22 

RM 

1 

0.232472 

0.85 

(1.1.2) 

1 

0.232489 

0.01 

(3,3,4) 

1 

0.232489 

(5,5,6) 

Klrchhoff 

0.232392 

2.04 

K 

3  dim 

0.232489 

We  see  that  for  hard  support  the  performance  of  the  model  (1,1,2),  X  =  1  is 
better  than  of  RM  model  with  X  -  0.87,  which  is  in  contrast  with  the  soft 


support.  Let  us  also  underline  that  the  energy  of  RM  model  with  X  =  0.87  is 
larger  for  hard  support  than  that  of  the  energy  of  the  3-dimensional  solution. 

5 

The  same  occurs  for  X  *  =•  and  soft  support. 

b 

In  Table  4.6  we  also  report  the  values  of  the  (nj.ng.i^)  model  when 
JfE 

a44  *  *  1~+  v  X  -  1  and  X  =  0.87.  We  see  that  it  is  advantageous 

for  the  model  (1,1,2)  but  detrimental  for  higher  models  compared  with  the 
solutions  of  the  3-dimensional  plate  model. 

In  Table  4.8  we  show  the  value  of  the  shear  force  on  the  diagonal 

Xj  =  Xg  and  the  error  when  compared  with  the  solution  of  the  3-dimensional 
plate  model. 


Table  4.8  The  shear  force  and  its  error  on  the  line  Xj 


x2  for  the 


soft  support. 


Xl  =  X2 

3  dim 

X  = 

=  1 

X  = 

0.91 

X  = 

0.87 

x  = 

=  5/6 

Q31 

€*/. 

Q31 

eV, 

Q31 

e*/. 

Q31 

e% 

0.4905 

-0.29 

-0.22 

24.9 

-0.26 

11.9 

-0.28 

5.3 

-0.29 

0.92 

0.4950 

-1.55 

-1.46 

6.0 

-1.  52 

2.3 

-1.55 

0.7 

-1.57 

0.77 

0.4968 

-2.78 

-2.79 

0.4 

-2.80 

0.8 

-2.81 

0.9 

-2.81 

0.81 

0.4986 

-4.75 

-5.01 

5.5 

-4.87 

2.6 

-4.  81 

1.3 

-4.75 

0.07 

0.49905 

-5.37 

-5.70 

6.  1 

-5.51 

2.5 

-5.41 

0.9 

-5.33 

0.76 

0.49995 

-6.70 

-7.08 

5.6 

-6.77 

0.9 

-6.62 

IB 

-6.49 

3.31 

0.5 

-6.79 

-7.  16 

5.5 

-6.  83 

0.7 

-6.  68 

■ 

-6.54 

3.53 

In  what  follows  we  always  use  in  the  RM  model  the  value  of  X  =  0.87, 
unless  otherwise  stated. 

In  Tables  4.9a  and  4.9b  we  show  the  moment  MJ2  computed  by  the  RM  (K  = 
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0.87)  and  (1,1,2)  (H  =  1.00)  model  for  the  hard  and  soft  support. 

Analogous  results  are  given  in  Table  4.2  for  the  three-dimensional  solution. 

In  the  table  we  also  show  that  relative  error  e  in  comparison  with  the 
three-dimensional  solution  whenever  it  is  larger  than  0.1%.  In  the  case  of 
hard  support,  both  RM  and  (1,1,2)  models  have  error  <  0.1%.  It  is 
interesting  that  the  model  (1,1,2)  produces  a  larger  error  for  the  soft 
support  moment  M12  than  the  RM  model.  Model  (3,3,4)  would  give  better 
results  than  the  RM  model. 


Table  4.9a  The  moment  ,  x^  =0.4  for  hard  and  soft  support. 


X2 

Xj  =0.4 

HARD 

SOFT 

RM 

(1,1,2) 

RM 

Crm% 

(1,1,2) 

Cl  12% 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.02368 

0.0018 

0.0018 

0.0018 

0. 

0.0018 

0. 

0. 11842 

0.0092 

0.0092 

0.0092 

0. 

0.0092 

0. 

0.21316 

0.0162 

0.0162 

0.0163 

0. 

0.0163 

0. 

0.35526 

0.0252 

0.0252 

0.0254 

0. 

0.0254 

0. 

0.45000 

0.0289 

0.0289 

0.0292 

0. 

0.0292 

0. 

0.48079 

0.0294 

0.0294 

0.0297 

0. 

0.0297 

0. 

0.49026 

0.0295 

0.0295 

0.  0286 

E23 

0.0283 

1. 1 

0.49500 

0.  0295 

0.0295 

0.  0239 

mm 

0.0246 

3.4 

0.49713 

0.0295 

0.0295 

0.0180 

0.0188 

5.03 

0.49903 

0.0295 

0.0295 

0.0087 

0.0085 

6.3 

0.49950 

0.0295 

0. 0295 

0.0045 

0.0047 

6.0 

0.50000 

0.0295 

0.0295 

0. 

Hi 

0. 

0. 
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Table  4.9b  The  moment  ,  x^ 


0.5  for  hard  and  soft  support. 


X2 

*1  ■ 

0.5 

HARD 

SOFT  | 

RM 

(1,1,2) 

RM 

(1,1,2) 

0. 

0. 

0. 

0. 

0. 

0.02368 

0.0019 

0.0019 

0. 

0.  | 

0. 11842 

0.0097 

0.0097 

0. 

0. 

0.21316 

0.0172 

0.0172 

0. 

0. 

0.35526 

0.0270 

0.0270 

0. 

0. 

0. 45000 

0.0315 

0.0315 

0. 

0. 

0.48079 

0.0323 

0.0323 

0. 

0. 

0.49026 

0.0324 

0.0324 

0. 

0. 

0.49500 

0.0325 

0.0325 

0. 

0. 

0.49713 

0.0325 

0.  0325 

0. 

0. 

0.49903 

0.0325 

0.0325 

0. 

0. 

0.49950 

0.0325 

0.0325 

0. 

0. 

0. 50000 

0.0325 

0. 0325 

0. 

0. 

In  Tables  4. 10a  and  4. 10b  we  show  analogous  data  for  the  shear  forces. 

Once  more  we  see  that  for  hard  support  the  accuracy  is  high  for  RM  and 
(1,1,2)  model.  Otherwise  we  see  once  more  that  the  values  are  better  for  RM 
model  (ft  =  0.87)  than  for  the  (1,1,2)  (ft  *  1)  model. 

We  will  see  in  Section  5  that  this  conclusion  is  not  necessarily  correct 
for  other  boundary  conditions. 
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Table  4.10a  The  shear  forces  Q^CXj.Xg) 


x 


*1  a 

ISM 

Q31 

Q32 

*2 

SOFT 

SOFT 

RM 

CRM*/- 

(1,1.2) 

Cn2% 

RM 

Crm5* 

(1,1,2) 

g112% 

0. 

0.419 

0. 

0.419 

0. 

0. 

■ 

0. 

— 

0.02368 

0.421 

0. 

0.421 

0. 

-0.638 

-0.685 

6.2 

0. 11842 

0.408 

0. 

0.408 

0. 

-3. 173 

-3.399 

6.2 

0.21316 

0.374 

0. 

0.374 

0. 

-5.617 

0.8 

-6.018 

6.2 

0.35526 

0.273 

0. 

0.273 

0. 

-8.844 

0.8 

-9. 474 

6.2 

0.45000 

0.  133 

0. 

0.  135 

1.5 

-10.363 

0.8 

-11.095 

6.2 

0.48079 

-0. 125 

— 

-0. 105 

— 

-10.559 

0.6 

-11.332 

6.5 

0.49026 

-0.807 

0.6 

-0.756 

6.9 

-9.844 

0.6 

-10.651 

7.5 

0. 49500 

-2. 188 

0.6 

-2. 157 

0.8 

-8.505 

1.  1 

-9.251 

7.5 

0.49713 

-3.573 

1.3 

-3.645 

3.4 

-7.569 

mm 

-8.221 

6.6 

0.49903 

-5.499 

0.6 

-5.815 

6.3 

-6.881 

HI 

-7.368 

6.0 

0.49950 

-6.073 

0.7 

-6.485 

6.5 

-6.721 

1.7 

-7.252 

6.4 

0.50000 

-6.678 

HI 

-7.221 

6.5 

-6.678 

IB 

-7.221 

6.5 

Table  4.11  The  value  of  Q^(0.5,  0.5)  for  soft  support  and  various  plate 
models. 


K 

Q31 

e% 

Model 

1 

-7.  15 

5.6 

RM 

0.91 

-6.83 

0.8 

RM 

0.87 

-6.68 

1.5 

RM 

0.845333 

-6.61 

2.8 

RM 

5/6 

-6.54 

3.4 

RM 

1 

-7.22 

6.5 

(1,1,2) 

1 

-6.82 

0.6 

(3,3,4) 

1 

-6.80 

0.3 

(5,5,6) 

0.87 

-6.73 

0.9 

(1,1,2) 

0.87 

-6.36 

6.3 

(3,3,4) 

0.87 

-6.34 

6.6 

(5,5,6) 

3  dim 

-6.97 

— 

We  see  that  for  the  RM  model  Jf  >  0.87  is  optimal.  For  the  (n,n,n+l)  model 
optimal  H  >  0.87  for  all  n  and  as  n  increases  the  optimal  H 
approaches  1  and  the  model  (1,1,2)  (K  =  0.87)  leads  to  smaller  error  than 
RM  model  with  H  =  0.87. 

In  Table  4.12  we  show  the  function  0^  for  the  RM  model  (H  =  0.87) 
and  (1,1,2)  model  with  K  =  1  which  are  defined  in  the  same  way  as  in  Table 
4.4  for  the  (three-dimensional)  plate  model. 
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Table  4.12  The  values  of  (S^tx^Xg),  x^  *  0.4,  0.5  for  various  models 
(soft  support) 


X2 

*1  s 

=  0.4 

*1  = 

mm 

RM 

(1,1,2) 

RM 

(1,1,2) 

0.49263 

3.25 

2.48 

2.21 

2.37 

0.49500 

3.24 

3.47 

2.23 

2.41 

0.49666 

3.24 

3.47 

2.20 

2.39 

0.49878 

3.23 

3.47 

2.04 

2.24 

0.49950 

3.23 

P.  47 

1.91 

2.  12 

0. 49989 

3.24 

3.47 

1.83 

2.03 

0.49991 

3.24 

3.47 

1.83 

2.  J3 

0.49994 

3.23 

3.47 

1.83 

2.03 

It  has  been  shown  in  [15]  and  [16]  that  when  the  boundary  is  smooth  the 
strength  of  the  boundary  layer  for  the  RM  model  is  P31  *  J  12H  =  3.23. 

The  detailed  theoretical  analyses  of  the  boundary  layer  behavioi  of  the 
(1,1,2)  model  is  not  available.  In  [17]  was  suggested  that  ^^(0.5,  x^)  ~ 

V  6X  =  2.28  for  the  RM  model. 

To  the  authors’  knowledge  there  is  no  rigorous  analysis  of  the  boundary 
layer  in  the  neighborhood  of  the  corners  of  u  for  RM  or  any  other  model. 

We  have  seen  that  as  d — >0  the  difference  between  the  soft  and  hard 
support  disappears  when  measured  in  the  energy  norm.  This  is  not  true  for 
the  values  Influenced  by  the  boundary  layer.  In  Table  4. 13  we  show  for  the 
RM  model  with  H  =  |  the  values  0^(0. 5,  0.5)  and  dQ^tO.S,  0.5)  for  d 
=  0.025  and  d  =  0.01  for  the  soft  support. 
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Table  4.13  The  values  of  0^(0. 5,  0.5)  sis  function  of  d 


d 

Q31(0.5,  0.5) 

dQ3i (0.  5,0.5) 

0.025 

-2.68 

-0.0671 

0.01 

-6.54 

-0.0614 

We  see  that  dQ^  converges  as  d — >0  as  expected  (but  not 
theoretically  proven).  We  can  expect  (see  [17])  that  for  the  RM  model 
CJ^^(0.5,  0.5)  is  proportional  to  J  it'.  In  Table  4.14  we  show  that  the 
expectation  is  correct. 

Table  4.14  The  dependence  of  0^(0. 5,  0.5)  on  H  for  the  RM  model. 

Q31(0.S,  0.5) 

-7.  15 
-6.83 
-6.68 
-6.54 

4.5  The  problem  of  the  L-shaped  domain. 

In  the  previous  Section  we  discussed  the  problem  of  the  square  plate. 
In  this  section  we  will  briefly  discuss  the  problem  of  the  simply  soft 
supported  L-shaped  plate.  The  domain  w  is  shown  in  Figure  4.7. 


0.  87 
5/6 


-7.  16 
-7.  16 
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Figure  4.7  The  L-shaped  domain. 

Once  more  we  will  consider  the  case  d  =  0,01  and  v  =  0.3.  We  will 
concentrate  here  on  the  behavior  of  the  solution  in  the  neighborhood  of  the 
origin  where  the  solution  is  singular.  We  will  present  only  the  results  of 
RM  model  for  H  -  0.87.  Table  4.5  shows  the  theoretical  strength  of  the 
singularities  for  the  RM  and  K  model  for  soft  support  (a  *  270°). 


Figure  4.8  The  shear  force  0^(0,  *2),  Q^tO,  x2^  for  the  ^  m°del. 
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Figure  4.8  shows  in  log  log  scale  the  shear  forces  on  the  line  x^  =  0  (as  a 
function  of  Xg).  In  the  figure  we  also  show  the  expected  theoretical  slope 
(0.3333). 

Figure  4.9  shows  the  values  of  the  moments  and  M12  as  well  as 

the  theoretical  slope.  We  see  that  the  singularity  is  different  for  the 
moments  and  the  shear  forces. 


iff7  lO"6  10"5  Iff4 


Figure  4.9  The  moment  MggCO,  x,,),  M12(0,  x2^  for  the  ^  model- 


11  we  show  the  shear  force  and  the  moment 

_2 

see  that  for  x2  >  10  (l.e.  of  the  thickness 

and  0.666  for  Mgg.  This  behavior  is  related 
to  the  Kirchhoff  singularity  (also  listed  in  Figure  4.5).  We  see  typically 
that  inside  the  domain,  in  the  distance  of  order  thickness  the  Kirchhoff 
model  describes  the  character  of  the  solution  very  well.  This  is  typical 
behavior  of  the  solution  of  every  plate  model  as  well  as  for  the  plate 
(three-dimensional)  problem. 


In  Figures  4. 10  and  4. 
in  a  larger  scale.  We 
the  slope  is  1.666  for 


) 
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SHEAR  FORCE  Q32  (0,  x2) 


5X104 

3X104 

1Y104 

5X103 

3X103 

1X103 
5x10  2 
3X10  2 

1X102 

5X10 

3xio 

lx  10 

5 

3 


10*6  lO"5  10"4  10"3  10‘2  10 


-1 


Figure  4.10  The  shear  force  Q^tO,  x2^ 


Figure  4.11  The  moment  M^O,  x2^ 
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Figures  4.10  and  4.11  show  the  behavior  inside  «.  Figures  4.12  and  4.13 
show  the  shear  force  QggtXj.  0)  (♦*»  reaction)  and  the  moment  M^x^  0). 


Figure  4.12  The  shear  force  QggCXj,  0) 


10 '7  10  -6  10  -5  10 A  10  -3  10  '2  10  -1 


Figure  4.13  The  moment  MjjtXj,  0). 
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We  see  that  the  singularity  in  the  close  neighborhood  of  the  vertex  is  the 
theoretical  one  as  follows  from  Table  4.5.  Nevertheless  the  Kirchhoff 
behavior  is  not  visible  because  the  line  Xg  =  0  is  on  the  boundary  and  the 
boundary  layer  is  not  captured  by  Kirchhoff  model.  The  shear  force 
O^Cx^O)  does  not  show  the  singularity  because  of  the  influence  of  the 
boundary  condition. 


Figure  4.14  The  shear  force  Q^Cx^O). 


We  have  shown  the  results  auid  asymptotic  behavior  for  the  RM  model 
only.  They  characterize  well  the  behavior  for  any  n-model  as  well  as  the 

4 

solution  of  the  three-dimensional  solution. 

We  have  seen  that  in  the  neighborhood  of  the  vertex  of  w  we  have 
exactly  two  kinds  of  singularities:  one  in  the  neighborhood  which  is  smaller 
than  the  thickness  and  auiother  one  —  the  Kirchhoff  —  which  is  valid  when 
the  distance  is  of  order  of  few  thicknesses. 
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In  Section  4.3  we  introduced  the  (classical)  Kirchhoff  plate  model  with 
B  k 

the  solution  w  =  Ug  .  The  Kirchhoff  model  is  not  able  to  distinguish 
between  the  hard  and  soft  simple  support.  Nevertheless  the  plate  problem 
(three-dimensional  as  well  as  other  plate  models)  converges  (in  the  scaled 
energy  norm)  to  the  Kirchhoff  solution  as  d — >0.  See  e.g.  [2],  [6]  and 
Theorem  4.  5. 

The  solution  u  of  the  Kirchhoff  model  satisfies  the  biharmonic 
equation 


n.2 

DA  w  =  q 


12(1  -  v) 


If  u>  is  a  polygon  and  simply  support  is  considered  the  problem  could  be 
also  formulated  as  follows:  Find  w  6  H2^)  such  that  (4.11)  holds  and 

(4. 12)  w  =  Aw  =  0  on  T  (except  vertices  of  w) 

This  formulation  is  equivalent  to  the  one  based  on  the  minimization  of 

g 

G(u).  Hence  we  can  formulate  the  problem  as  the  system 

Aw  -  v 

(4. 12)  Av  =  ^ 

w  *  v  =  0  on  T  (except  on  the  vertices  of  w) 

and  w  c  H^tu).  The  conditions  that  w  e  u)  is  essential.  When  w  is  a 
convex  polygon  then  u  can  be  found  directly  by  (4.12).  In  this  case  v  e 
H^u)  and  w  €  H 2(w).  Nevertheless,  when  a  >  n  then  the  solution  of 
(4.12)  yields  v  €  H^u)  and  w  €  H^w)  but  w  t  \?{w)  and  so  we  have  to 
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add  (respectively  subtract  a  singular  function  e  Hitt))  so  that  w  e 
H 2(o»).  We  have  now 


For  a  <  n,  a  *  g  : 


w  -  Cr  0(0)  +  smoother  terms 


For  a  >  n  : 


when  +  2  <  —  ,  then 
^  a 


w  =  Cr  7r/a  +  20(0)  +  smoother  terms 


when  —  <  --  +  2  ,  then 
a  a 


w  *  Cr2TL/a0(0)  +  smoother  terms 


when  —  *  --  +  2  ,  then 
a  a 


w  *  r2lI/a[ ic^ce)  +  c2^2(0) 


Here  r,0  are  polar  coordinates  with  the  origin  in  the  vertex  of  the  domain 

w  and  ^(9)  is  a  smooth  function.  The  case  a  =  is  an  exceptional  angle 
5-e 

and  then  w  e  H  (w),  e  >  0,  arbitrary. 

The  moments  M^,  i,J  =  1,2,  and  the  shear  forces  Q^,  i  =  1,2,  are 
given  in  the  standard  form 
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+  V 


(4.13) 


M 


11 


”22 


2 

aw 


+ 


) 

) 


M  =  D(l-y)[  — —  ] 

l  axx  ax2  J 


Si  =D(-^+  ~~2  ] 

*■  3xi  axtSxi  j 


°02  ■  D( 


a3w  +  a' 


2 

axi  ax2 


a3w  "I 

ax3  J 


and  hence 


M,  *  Cr*1  ,(0)  +  smoother  terms 


U 


i.J 


Q^j  -  Cr*2  Vj(0)  +  smoother  terms 


As  before  the  coefficients  and  A 2  for  the  Kirchhoff  model  are  given  in 

Table  4.5. 

As  was  previously  mentioned,  the  Kirchhoff  model  cannot  distinguish 
between  hard  and  soft  support.  We  show  now  that  the  Kirchhoff  model  should 
be  understood  as  an  approximation  of  the  hard  support  although  sometimes  it 
is  used  for  the  soft  support  also  (see  below). 

In  Table  4. 15  we  show  the  error  of  the  Kirchhoff  model  when  compared 
with  the  solutions  of  the  (three-dimensional)  plate  problem  of  the  square 
uniformly  loaded  plate  (v  =  0.3,  d  =  0.01)  for  d  =  0. 1  and  d  =  0.01. 
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Table  4.15  Comparison  of  the  energy  of  the  Kirchhoff  model. 


d 

f,  B,  .  A,  H.  ,1 

|e  (w)  -  c  (u  ) | 

1/2 

|eB(w)  -  cA(uS) | 

1/2 

A ,  S. 
e  (u  ) 

W  J 

A,  S, 
e  (u  ) 

Em 

20.31*/. 

39.56% 

0.01 

2.03% 

11.87% 

In  Table  4.16  we  show  for  d  =  0.01  some  comparison  of  data  from  the 
Kirchhoff  model  and  the  three-dimensional  solution. 


Table  4. 16  Comparison  of  shear  forces  and  moments  of  the  Kirchhoff  and 


three-dimensional  model  (d  =  0.01). 


Q3i(0.4,  0.5) 

Q3l(0.5-0,  0.5) 

M12(0.5-0,  0.5) 

hard  support 

0. 

0. 

0.0325 

soft  support 

-9.72 

-6.79 

0. 

Kirchhoff 

0. 

0. 

0.0325 

In  the  table  we  denoted 


,(0.5-0,  0.5)  =  lim  Q(x1t  0.5) 
0,1  Xi<0.5  1 

xi  ->0.5 


Kirchhoff  model  also  reproduces  the  paradox  of  the  hard  simple  support  i.e 

we  have  w  — >w  as  m — >»  (see  Theorem  4.2)  and  [6],  [13]  and  [14]. 

m  c 

Often  in  the  Kirchhoff  theory  the  reaction  e.g.  V2  =  (for  x2 

0.5)  is  computed  from  the  formula 

(414)  V2  ’  (^2  *  ife) 


46 


see  e.g.  [18].  This  derivation  tries  to  simulate  the  reaction  for  the  soft 
simple  support  when  using  essentially  the  data  from  the  hard  support.  The 
main  idea  of  the  standard  derivation  of  (4.14)  is  following:  The  virtual 
work  by  the  Kirchhoff  theory  can  be  written  in  the  form 

(4.15)  B(w,v)  =  J  qvdx  -  <f  Qnvdt  +  +  Mnt^Jdt 

w 

Here  we  denoted  by  Mn  the  normal  moment  (for  the  simple  support  Mn  =  0) 
and  M  ^  is  the  twist  moment.  By  integration  by  parts  in  (4.15)  we  get 

(4.16)  B(w,v)  =  J  qvdx  -  <f  [q„  +  ^fjvdt  +  |  M^t 

w  r 

Here  is  (generalized)  derivative  with  respect  to  the  tangent.  We 

remark  that  Mnt  is  discontinuous  by  (in  the  case  of  the  square  domain) 

piecewise  smooth  and  hence  — is  a  Dirac  function  (concentrated  force)  in 

at 

the  vertex  with  the  value  of  the  "jump"  in  the  moment.  In  the  general  case, 
for  example  when  a  >  then  the  "Jump”  is  infinite  and  so  it  is  not 
completely  plausible  to  make  a  reasonable  physical  interpretation.  The  idea 
in  deriving  (4.16)  is  to  remove  the  virtual  work  by  the  twist  moment  but 
still  keep  the  same  value  of  the  total  virtual  work.  In  the  case  of  the 
square  plate  (a  *  90°)  the  concentrated  force  computed  through  the  twist 
moment  is  relatively  successful.  We  get  here  ft  =  0.0325  *  M12(0.5,  0.5) 
computed  with  the  exact  value  ft  defined  in  Section  4.2  (ft  *  0.0290)  and 
the  difference  is  about  10%. 

4.7  Summary  and  conclusions. 


The  simple  support  is  an  idealization  which  does  not  necessarily 


describe  the  reality  well.  There  are  many  possible  formulations.  We 

mentioned  only  two  possibilities,  the  hard  and  soft  simple  support.  There 

are  others.  For  example  we  mention  the  notion  of  the  supersoft  support  which 

-d/2 

is  defined  so  that  H(G)  =  {u  |  I  u^Xj.Xg.x^dXg  =  0,  (x^Xg)  e  D 

-d/2 

Xj.Xg.Xg)  =0  as  was  assumed.  This  will  Influence  the  solution 
by  still  various  basic  features  discussed  earlier  for  the  soft  support  will 
occur. 

In  the  literature  sometimes  another  formulation  of  the  soft  support  is 
used.  Here  X(Q)  =  {u  luCx^.Xg.O)  =  0,  (Xj.Xg)  €  T>.  Nevertheless  this 
formulation  has  no  sense  because  for  the  3-dimensional  formulation  the 
quadratic  functional  G(u)  will  be  not  bounded  from  below  and  the  solution 
does  not  exist.  (The  essential  reason  is  that  under  a  concentrated  load 
(reaction)  the  displacement  is  infinite. ) 

Comparing  the  results  we  mentioned  earlier  we  see: 

(1)  The  Kirchhoff  model  gives  completely  unreliable  results  in  the 
neighborhood  of  the  boundary  (of  size  of  the  plate  thickness) 
although  it  gives  acceptable  results  inside  the  domain.  It 
approximates  much  better  the  hard  support  than  the  soft  support. 

In  the  close  neighborhood  of  the  vertices  the  Kirchhoff  model  leads 
to  the  singularities  which  are  very  different  when  compared  with 
the  exact  solution  of  the  3-dlmensional  problem. 

(2)  There  is  a  significant  difference  between  the  solution  for  the  hard 
and  soft  simple  support.  This  difference  is  limited  to  the 
neighborhood  of  the  boundary.  The  size  of  the  neighborhood  depends 
on  the  smoothness  of  the  boundary.  The  soft  support  does  not  have 
certain  paradoxical  properties  which  the  hard  one  hats. 


instead  Ug( 
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(3)  The  RM  model  performs  well.  The  shear  factor  J <  has  a  positive 
influence  on  the  quality  of  the  solution.  Nevertheless,  the 
optimal  value  depends  on  the  concrete  setting  and  aims  of 
computation.  In  the  given  case  the  value  H  -  0.87  is  close  to 
optimal. 

(4)  The  singularity  of  the  (1,1,2)  model  is  the  same  as  of  the 
three-dimensional  solution.  The  RM  model  leads  to  a  slightly 
different  singularity. 

(5)  The  singular  behavior  of  the  RM  model  and  the  n-model  is  well 
described  by  the  theory  which  is  valid  in  the  area  of  approximately 
^  +  1  thickness.  In  the  areas  of  larger  distance  the  solution  has 
singular  behavior  well  described  by  the  Klrchhoff  model. 

(6)  The  RM  model  (X  *  0.87)  performs  better  than  the  model  (1,1,2) 
with  X  *  1  (not  with  X  *  0.87)  but  the  model  (3,3,4)  with  X  = 

1  outperforms  RM. 

(7)  The  negative  reaction  of  the  solution  of  the  three-dimensional 
solution,  RM  model,  and  n-model  differs  approximately  10%  from  the 
Kirchhoff  negative  concentrated  reaction. 

5.  The  plate  problem  for  various  boundary  conditions. 

In  Section  4  we  addressed  in  detail  the  problem  of  the  square  simply 
supported  plate.  In  this  section  we  will  briefly  analyze  the  essentials  when 
other  boundary  conditions  are  Imposed  especially  with  respect  to  the 
differences  to  the  case  discussed  earlier. 
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5. 1  The  square  clamped  plate. 


Let  us  consider  the  square  plate  with  all  build  in  boundary  conditions. 
Here  W(Q)  =  (u  |  u  =  0  on  S>.  As  before,  we  can  similarly  define  RM 
model,  n-model,  and  Kirchhoff  model. 

Analogous  to  our  previous  analysis,  we  use  the  matrix  B  and  parameter 
H  for  the  RM  model.  Then  as  d — >0  the  difference  (measured  in  the 
scaled  energy  norm)  converges  to  zero.  Hence  the  problem  has  quite  analogous 
properties  as  before. 

For  the  numerical  analysis  we  still  consider  as  before  the  unit  square 
plate,  d  =  v  =  0.3  and  for  the  RM  model  we  will  use  H  =  0.87,  for 

the  n-model  we  will  use  H  -  1.  In  Tables  S. 1  and  5.2  we  report  the  moment 


M11^X1,X2^  3111(1  M22^X1,X2^  for  X1  =  x2*  0  2777  respectively.  The 

moment  M^tx^O)  is  essentially  identical  (e  <  0.1%)  for  all  three  models 
considered.  Hence  no  error  is  present.  On  the  other  hand,  the  error  in 


M^CXi.O)  for  the  RM  model  is  large  (30%),  while  for  the  (1,1,2)  model  it  is 
acceptable.  Here  we  see  a  significant  difference  in  comparison  to  the  simple 


support.  Let  us  mention  that  the  Kirchhoff  plate  model  yields  M^CO,  0.5)  = 
0.0513,  Mpp,(-0,  0.5)  =  0.0153,  i.e.  the  identical  results  as  the  RM  model. 

Moment  M^  has  a  boundary  layer  which  is  not  captured  by  the  RM  model. 
As  in  Section  4.2  (see  Table  4.4)  we  define  ^xl,x2^  so  that 


|M22(Xi,  X2)~M22>(Xi,  X2)  | 

|M™(0.5,  X2)  -  M22>(0. 5,  x2)| 


where  the  function  022^xl,x2^  characterizes  the  boundary  layer.  We  use  in 
the  above  formula  M^tXj.Xg)  as  natural  smooth  extension  of  M^2^X1’X2^ 
from  the  inside  where  =  M^J.  In  Table  5.3  we  show  a  few  values  of 


&1-2'  -  4dl" 
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Table  5.1  The  moments  M^tXj.O),  for  different  plate  models. 


Mn(xi,0) 

M2z(Xi »  0) 

*1 

3  dim 

RM 

(1,1.2) 

3  dim 

RM 

CRM5* 

(1.1,2) 

Cii2y‘ 

0. 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

0. 

-0.0229 

0. 

0.2000 

-0.0157 

-0.0157 

-0.0157 

-0.0163 

-0.0163 

0. 

-0.0163 

0. 

0.4000 

0.0164 

0.0163 

0.0164 

0.0026 

0.0027 

0. 

0.0026 

0. 

0.4900 

0. 04r<  J 

0.0470 

0.0470 

0.0141 

0.0141 

0. 

0.0141 

0. 

0.4930 

0.0483 

0.0483 

0.0483 

0.0145 

0.0144 

0. 

0.0144 

0. 

0.4990 

0.0509 

0.0509 

0.0509 

0.0168 

0.0152 

9.2 

0.0170 

1.4 

0.4995 

0.0510 

0.0510 

0.0510 

0.0176 

0.0153 

12.3 

0.0179 

2.  1 

0.4999 

0.0512 

0.0512 

0.0512 

0.0207 

0.0153 

25.3 

0.0212 

2.8 

0.5000 

0.0513 

0.0513 

0.0513 

0.0220 

0.0153 

30.0 

0.0220 

0. 

models. 


MxiIxl  0.2777) 

M22(xi, 0.2777) 

*1 

3  dim 

RM 

(1,1,2) 

3  dim 

RM 

(1,1,2) 

e  112°^ 

0. 

-0.0102 

-0.0102 

-0.0102 

-0. 0075 

-0.0075 

0. 

-0.0075 

0. 

0.2000 

-0.0075 

-0.0075 

-0.0075 

-0.0058 

-0.0058 

0. 

-0.0058 

0. 

0.4000 

0.0083 

0.0083 

0.0083 

0.0015 

0.0015 

0. 

0.0015 

0. 

0.4900 

0.0256 

0.0256 

0.0256 

0.0077 

-0.0077 

0. 

0.0077 

0. 

0.4930 

0.0263 

0.0263 

0.0263 

0.0079 

0.0079 

0. 

0.0079 

0. 

0.4990 

0.0279 

0.0279 

0.0279 

0.0092 

0.0084 

9.2 

0.0093 

1.4 

0.4995 

0.0280 

0.0280 

0.0280 

0.0096 

0.0084 

9.0 

0.0098 

0.4999 

0.0281 

0.0281 

0.0281 

0.0113 

0. 0084 

25.3 

0.0116 

iiil 

0.5000 

0.0281 

0.0281 

0.0281 

0.0121 

0.0084 

30.  1 

0.0121 

0.  1 
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Table  5.3  The  coefficient  for  the  clamped  plates. 


X1 

3  dim 
022[xl • 0) 

(1,1,2) 

022(X1>O) 

0.4990 

15.43 

1 

0.4996 

18. 17 

0.4999 

23.85 

13.32 

0.49993 

24.  13 

13.31 

0.49999 

24.78 

13.29 

We  see  that  the  RM  model  does  not  have  essentially  any  boundary  layer  for  the 
clamped  boundary  condition.  This  follows  from  [15],  [16].  On  the  other  hand 
the  model  (1,1,2)  and  the  plate  problem  shows  a  strong  boundary  layer.  The 
theoretical  strength  of  the  solution  of  (1,1,2)  model  is  ~  13-1 

(see  [19])  which  is  very  close  to  the  values  shown  in  Table  5.3.  The  exact 
solution  shows  still  stronger  boundary  layer. 

5.2  The  plate  with  two  opposite  sides  clamped  and  two  free. 

in  this  section  we  will  consider  the  case  when  the  type  of  the  boundary 
condition  is  changing  in  the  vertex.  We  will  consider  the  unit  square  domain 
with  the  sides  x2  =  ±0.5  to  be  clamped  and  x1  =  ±0.5  to  be  free.  The 
boundary  condition  in  the  vertex  A  *  (0.5,  0.5)  changes  the  type  from  one 
side  to  the  other.  Once  more  as  d — >0  the  solution  of  the  plate  problem 
(three-dimensional)  as  well  as  the  n-model  solution  converges  in  (scaled) 
energy  norm  to  the  Klrchhoff  solution.  See  for  example  arguments  in  [2], 

[6].  In  the  neighborhood  of  the  vertex  A  the  solution  has  singular 
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behavior.  Denoting  by  r  the  distance  from  the  vertex  A  we  have 

*  Ci1'*1  Vjj(0)  +  smoother  terras 
Qgj  =  Cgr*2  Vj(®)  +  smoother  terms 

Where  A^  and  X^  are  different  for  different  models.  Their  values  for  our 
case  sire  given  in  Table  5.4. 

The  value  X2  =  1  for  the  RM  and  (1,1,2)  model  means  that  for  this 
singularity  the  angle  is  exceptional  and  the  shear  force  C^j  is  bounded 
in  the  neighborhood  of  A. 


Table  5.4  The  singularity  coefficients  Aj,  X 2  for  various  models. 


RM 

(1,1,2) 

3  dim 

K 

*1 

0.7583 

0.7112 

1.0687 

+ 

i  0.4386 

X2 

1. 

1. 

0.687 

+ 

i  0.4386 

Let  us  now  consider  the  case  as  in  the  previous  section,  namely  d  = 

on  the  line  x^  =  x 2  of  the 
(three-dimensional)  solution  of  the  plate  problem  sis  a  function  of  the 
distance  r  from  the  point  A. 


0.01,  v  =  0.3.  Figure  5.1  shows  moment  Mgg 
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three-dimensional  solution. 


In  Figure  5. 1  we  also  show  the  theoretical  singularity  -0.289  as  follows 

-2 

from  Table  5.4.  For  r  >  10  the  moment  shows  the  character  of  the 
Kirchhoff  singularity  where  the  singularity  coefficient  is  complex. 

In  Figure  5.2  we  show  in  more  detail  in  log  log  scale  the  moment  for  the 
RM  model  and  three-dimensional  solution.  We  clearly  see  different  growth  as 
follows  from  Table  5.4. 


Figure  5.3  shows  the  behavior  of  the  moments  Mjj,  and  and  shear 

force  Qgg  from  the  RM  model.  We  clearly  see  the  growth  in  accordance 


with  Table  5.4. 


As  we  have  seen  in  Figure  5.1,  the  moment  Mg2  shows  oscillations.  To 
understand  this  oscillation  better  we  show  in  Figure  5.4,  the  moment  of 

the  RM  model  as  a  function  of  Lg  r.  In  addition  we  show  the  function 
Ccos(0.4386  lg  — )  as  approximations  by  the  function  having  the  Kirchhoff 

ro 

type  singularity.  We  see  very  good  agreement. 
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MOMENT 


Figure  5. 


2  The  moment  ^22^*"*^  ^e  RM  model  and  three-dimensional 

solution.  . 


10 '5  10  A  10  *3 


30 

20 

10 
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Figure  5.3  The  moments  M 


IT  M 22  the  shear  force  Q32  of  the  m  model 


SHEAR  FORCE 


Figure  5.4  The  moment  of  RM  solution  and  its  approximation  by 

Kirchhoff  type  singularity. 


Table  5.5  The  moment 


*>2  °n  X1 


Xg  for  various  models. 


*1 

3  dim 

RM 

n  =  o.9i 

cV. 

RM 

H  =  0.91 

cV. 

RM 

H  =  5/6 

e*/. 

(1,1,2) 

e% 

0. 

-0.0407 

-0.0407 

0. 

-0.0407 

0. 

-0.0407 

0. 

-0.0407 

0. 

0.  135 

-0.0317 

-0.0317 

0. 

-0.0317 

0. 

-0.0317 

0. 

-0.0317 

0. 

0.27 

-0.0048 

-0.0048 

0. 

-0.0048 

0. 

-0.0048 

0. 

-0.0048 

0. 

0.36 

0.0239 

0.0238 

0. 

0.0238 

0. 

0.0238 

0. 

0.0238 

0. 

0.45 

0.0612 

0.0612 

0. 

0.0612 

0. 

0.0612 

0. 

0.0612 

0. 

0.4635 

0.0670 

0.0670 

0. 

0.0670 

0. 

0.0670 

0. 

0.0670 

0. 

0.4815 

0.0722 

0.0722 

0. 

0.0722 

0. 

0.0722 

0. 

0.0722 

0. 

0.4953 

0.0624 

0.0626 

0.2 

0.0627 

0.4 

0.0627 

0.4 

0.0621 

0.5 

0.4968 

0.0575 

0.0576 

0.  1 

0.0578 

0.5 

0.0579 

0.7 

0.0567 

1.3 

0.4995 

0.0630 

0.0624 

0.9 

0. 0632 

0.4 

0.0640 

1.6 

0.0605 

4.0 

0.49995 

0. 1173 

0.  1125 

4.0 

0.  1140 

2.7 

0.  1154 

1.6 

0. 1143 

2.5 

0.499995 

0.2252 

0.  1989 

11.7 

0.2014 

10.5 

0.2038 

9.5 

0.2212 

1.7 
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In  Table  5.5  we  show  some  values  of  the  moment  Mgg  for  various  models. 
We  see  that  for  x  *  0.5  the  moments  of  the  model  (1,1,2)  are  better  than 
for  the  RM  model.  It  relates  directly  to  the  different  singular  behavior  of 
the  solution  from  the  various  models,  as  follows  from  Table  5.4. 

5.3  Summary  and  conclusions. 

Let  us  mention  some  conclusions  which  follow  from  cur  analysis: 

1)  The  boundary  layer  behavior  of  the  RM  model  as  described  in  [15] 
[16]  is  very  weak  for  the  clamped  boundary  condition.  Nevertheless 
the  strong  boundary  layer  is  present  for  the  three-dimensional 
formulation  as  well  as  (1,1,2)  model.  The  RM  model  leads  to  a  very 
large  error  here  when  compared  with  the  exact  solution  of  the 
three-dimensional  problem.  In  contrast  the  model  (1,1,2)  gives 
very  acceptable  results. 

2)  Changing  the  type  of  boundary  conditions  in  the  vertex  usually 
leads  to  different  singular  behavior  of  the  solution  of  the  RM 
model  and  three-dimensional  solution.  The  model  (1,1,2)  has  the 
same  strength  of  the  singularity  as  the  three-dimensional  solution. 
Hence  the  RM  model  yields  a  large  error  in  a  very  small 
neighborhood  of  the  vertex.  The  structure  of  the  solution  is 
further  complicated  by  the  different  strengths  of  boundary  layer  on 
the  sides  with  different  boundary  conditions. 

6.  Additional  remarks  and  conclusions. 

The  solution  of  the  plate  problem  is  very  sensitive  to  boundary 
conditions  in  the  neighborhood  of  the  boundary.  Different  models  lead  in 
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general  to  very  different  boundary  layer  behavior  and  singularities  in  the 
neighborhood  of  the  vertices  of  the  domain.  We  discussed  only  the  behavior 
of  the  moments  and  shear  forces,  i.e.  of  the  resultants,  the  detailed 
solution  (stresses)  has  still  more  complicated  structure  (see  e.g. ,  [20]). 

The  singular  and  boundary  layer  behavior  of  course  strongly  influences 
the  accuracy  of  the  numerical  solution  and  requires  very  careful  mesh  design. 

If  the  plate  is  relatively  thin,  say  of  the  diameter,  the 

sensitivity  of  the  data  inside  of  w  is  much  smaller  than  in  the  boundary 
layer.  Nevertheless  it  still  can  be  large. 

The  design  of  the  model  has  to  be  directly  related  to  the  alms  of 
computations.  For  example  the  Kirchhoff  model  typically  could  lead  to 
results  which  are  very  far  from  the  exact  three-dimensional  solution  if  the 
data  at  the  boundary  are  of  interest.  (These  data  are  reported  in  any  book 
about  plate  theory.)  Inside  of  the  domain  the  Kirchhoff  data  usually  are 
relative  reliable  but  it  depends  on  the  type  of  boundary  conditions  and  the 
structure  of  the  plate  domain. 

Reissner-Mindlin  model  performs  well  and  captures  well  the  boundary 
layer  behavior  for  some  boundary  conditions;  nevertheless  it  fails  completely 
for  some  others  such  as  clamped  boundary  conditions  and  could  yield  results 
with  error  of  3054.  The  (1,1,2)  model  especially  with  a  shear  factor  is  much 
more  relevant  and  gives  reliable  data  also  when  RM  fails  as  in  the  case  of 
the  clamped  boundary  condition.  This  model  also  has  the  same  singular 
behavior  as  the  solution  of  the  (three-dimensional)  plate  model,  while  the  RM 
often  has  different  singularity.  The  singularity  of  the  solution  of  the 
Kirchhoff  model  is  completely  different  when  compared  with  the  singular 
behavior  of  the  solution  of  three-dimensional  problem. 

The  selection  of  the  plate  model  has  to  be  related  to  the  aims  of 
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computation  and  the  best  way  is  to  select  it  in  an  adaptive  feature  from  a 
family  of  hierarchical  models;  in  different  plate  areas.  In  addition 
a-posteriori  error  of  the  solution  when  compared  with  the  exact  three- 
dimensional  solution  is  desirable  and  can  be  made  by  comparing  the  results 
from  different  hierarchically  ordered  models. 

Finally  we  mention  that  the  boundary  layer  influences  the  accuracy  of 
the  finite  element  solution  and  has  to  be  dealt  with.  (For  example  by  an 
adaptive  solver.  Furthermore,  the  finite  element  solution  of  the  n-model 
leads  to  the  locking  effects  which  are  dealt  with  in  various  ways.  In  [21] 
we  have  shown  that  the  p-version  of  FEM  for  p  k  4  does  not  show  practically 
the  locking  effects. 
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